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ABSTRACT 

The purpose of this report is to analyze the load carrying capacity of an 
external pressurized gas journal bearing including the hydrodynamic effect 
caused by the rotational speed of the journal. The analysis is an approximate 
first order pertubation solution, i.e. it assumes small bearing eccentricity 
ratio. Furthermore, it is assumed that the gas feeding takes place through 
orifices in the centerplane of the bearing and that the number of feeding holes 
is sufficiently large to be considered a line source. 

Results are given for the load carrying capacity and for the attitude 
angle. In addition an investigation is made of the validity of the appraxima- 


tion used in arriving at the solution. 


INTRODUCTION 

Gas journal bearings are principally of twotypes: 1) self acting, and 
2) externally pressurized. The self acting bearing works by hydrodynamic 
action and, therefore, has limited load capacity. It requires a small 
clearance and tight manufacturing tolerances. The externally pressurized 
bearing relies on hydrostatic pressure to generate its load carrying 
capacity. It is used in applications with higher load requirements and with 
less manufacturing precision but requires an external pressure supply. 

Although the externally pressurized bearing is primarily hydrostatic, 
the rotation of the journal causes some hydrodynamic action. This effect is 
in most cases negligible because of the relatively large clearance but may 
become significant at higher speeds and at large eccentricity ratio. Since 
in recent times such cases have been encountered it is of value to develop 
an analysis for the combined hydrostatic-hydrodynamic gas journal bearing. 


This is the purpose of the present report. 


RESULTS 
The results of the analysis are presented in Figures l-l2 where curves 
for load carrying capacity and attitude angle are given. Four dimensionless 


numbers are necessary for a bearing calculation: 
the bearing number A = hue (BY 
a 


the bearing length-to-diameter ratio = L/D 


the feeding parameter Ay = Sadie [Bx 


the pressure ratio = Psupply (Pent 


An estimate of the validity of the results can be obtained from Figures 
21-28 where the dotted curves are the previous results and the solid curves 
present an improved analysis. The improved analysis, however, only 
satisfies the analysis within a certain error which is marked on the curves. 
Thus a comparison between the two analyses is only valid when the error 


is small. 


spoereny 
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DISCUSSION OF ANALYSIS 

Theoretical analyses exist separately for the self acting bearing and for 
the externally pressurized bearing under the assumption that the eccentricity 
ratio@ is small (Ref. 2 and 3). Both these solutions are first order pertuba- 
tion solutions around € = 0 such that the bearing load becomes proportional 
to €. The results are in good agreement with more extensive solutions and 
with experimental results. Thus it seems logical to attempt an analysis that 
combines the two solutions. 

Assuming isothermal conditions in the gas film the exact equation for the 


problem under consideration is Reynolds' equation: 


SR [oh HE] + RL 5H BE] = bu Reo Y - 2pm 


where P = film pressure, 9 = gas density, h = film thickness, R = journal 
radius, W = angular speed of journal, ss = gas viscosity, x = circumferen- 
tial coordinate, z = axial coordinate and m = mass flow per square inch from 
the external pressure source (m is a function of the local pressure P). This 
mass flow enters the equation only at the discreet points at which flow is 
introduced. For this reason it is extremely difficult to obtain a solution of 
the above equation even in approximate form. To circumvent this difficulty 
an assumption shall be made in regard to the nature of the gas feeding. The 
gas feeding takes place through orifice restricted feeding holes in the center- 
plane of the bearing. It shall be assumed that the number of feeding holes is 


sufficiently large to be treated as a line source. Thus the gas feeding no 
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longer enters the differential equation but instead becomes a boundary condi- 
tion.to the equation. This boundary condition stipulates that at the centerplane 
of the bearing the pressure gradient must be such that the flow into the 
bearing is equal to the outflow from the feeding holes. The flow through the 
feeding holes is restricted by orifices. Any restriction due to the annular 
flow area existing between the edge of the feeding hole and the journal surface 
is neglected. 

The normal equation for flow through an orifice has a discontinuity when 
the pressure ratio reaches a critical value. This discontinuity is eliminated 
by using a modified orifice flow equation which has proven to give satisfactory 
results. 

Based on the above assumptions the equation defining the gas feeding can 


be set up. It is found to be governed by the dimensionless feeding parameter: 
= SaNoa t/ Rr. 
A, ~ RC m 


(m= viscosity, P, = ambient pressure, C = radial clearance, N = number of 
feeding holes, & = orifice coefficient, a = orifice radius, m = number of 


orifices in series, R= gas constant, T = total temperature. ) 


and the ratio of supply pressure to ambient pressure: 


Returning to Reynolds' equation it is written in dimensionless form disclosing 
two additional dimensionless parameters: the bearing number 
2 
A= Sae(B) 
R \C 
(w = angular journal speed, R = journal radius) 
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and the bearing length-to-diameter ratio: 

B= 5 
The solution will be a function of these 4 parameters. A first order pertuba- 
tion solution of Reynolds' equation is desired. 

Thus the film pressure is approximated by: 

Pe R+eP 
P,, is the pressure distribution when the journal] is concentric in the bearing, 
and is found to be: 

B= Vi+ gxe-t) 
where t is a constant determined by the bearing-feeding system and (E-T) 
is the axial distance measured from the end of the bearing, i.e. pe is a 
linear function. Substituting the above first order expression for the pressure 
into Reynolds' equation yields the first order pertubation equation with 
boundary conditions, see Appendix, eq. 5: 

Hho) + dBp) — 4 MEE) = AR sing 

Unfortunately this equation has no simple solution due to the form of P,: 


Instead the nature of the solution is explored. This is done by approximating 


ee 


where P, ae is the average value of Po: Then the equation is similar to 


+ by a constant: 


Ref. 1 and 2 and a solution can be obtained. Numerical values of the dimen- 
sionless load carrying capacity and the attitude angle are given in Figs. 1-12 


as a function of A for a wide range of A,-values, two L/Deratios and 
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three pressure ratios. 

The second part of the analysis evaluates the validity of the approxi- 
mation + = constant and establishes the regions in which the results are 
valid. Whereas previously + was constant throughout the bearing, it 
now varies along the length of the bearing. Since Pp is linear it is 
reasonable to replace t with a linear or almost linear expression. 
This leads to a solution involving Bessel functions of a complex variable 
which is rather impractical to evaluate and instead the Bessel functions 
are approximated by complex trigonometric functions. Hence the solution 
does not satisfy the differential equation exactly and the discrepancy is 
presented in Figures 13-20. This error must be small compared to 1. 
Otherwise the analysis proceeds as before and results for load carrying 
capacity and attitude angle are given in Figures 21-28. Inthe same figures 
are also shown the results for 4 = constant for comparison. In addition 
the error in using the trigonometric approximation is indicated on the 
curves as taken from Figures 13-20. For A less than 1.1 the agreement is 
good for the load carrying capacity and reasonably good for the attitude 
angle. No obvious correlation can be found between the magnitude of the 
calculated error and the point of divergence of the results. For the parti- 
cular case of L/D = 2, V =2, and A, = .5, the error never exceeds 2% 
and even so the results of the two analyses differ considerably. This points 
to the need for future investigation. 


Finally, it shall be investigated if the total load carrying capacity can be 
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considered as a simple superposition of the hydrodynamic and the hydre- 
static load. Figures 2-24 indicates that this is not the case. However, 
since Af Sayer (BY , where P, can be interpreted as being the mean 
film pressure for a purely hydrodynamic peeing. AL should be replaced 
by A & in applying the superposition for the hydrodynamic-hydro- 
static bearing. This has been done in Figures 29-32. It is seen that even 
with the modified A the superposition still is not valid. This is not 
surprising since the hydrodynamic and the hydrostatic action are coupled 


in a "non-linear" way through the feeding conditions. 


CONCLUSIONS 
The load and the attitude angle of the hybrid bearing are determined by 
4 dimensionless numbers, A , L/D, Ay and P,/P,. 
The rotation of the shaft increases the load carrying capacity and the 
attitude angle of the externally pressurized bearing. The effect is 
most pronounced for high values of the bearing number A (i.e. at high 
speeds), for low values of the supply pressure and for high values of 
the bearing length-to-diameter ratio. 
The load carrying capacity of the hybrid bearing is not a superposition 
of the hydrostatic and the hydrodynamic load even if the bearing number 
A. is based on the average film pressure instead of, as usual, the 


ambient pressure. 


RECOMMENDATIONS 
The analysis should be extended to obtain the complete solution of 
the first order pertubation equation. 
An investigation of the effect of higher values of the eccentricity ratio 
should be undertaken. 
The stability of the hybrid bearing should be studied especially to find 
out how hydrostatic pressure influences the stability of hydrodynamic 


bearings. 
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NOMENCLATURE 


Constants specifically defined in the analysis for the purpose of sim- 
plification only are not listed below. 


Orifice radius inch 
Cc Radial bearing clearance inch 
Cy d Constants given by equation (20) - (21) 
x (1st part) and (55) - (54) (2nd part) 
F F ; Radial and tangential force components lbs. 
G Complex form of H, see equation (12a) 
s,(9), s,($) Functions defined by equation (15) - (16) 
( 1st part) and (48) - (49) (2nd part) 
H eP P 
ae | 
h Film thickness 2, inch 
h Dimensionless film thickness = = 
L Bearing length inch 
M Mass flow through one feedinghole lbs. sec 
in 
m Number of orifices in series in one 
feeding hole 
N Number of feeding. holes around circum- 
ference 
P Pressure in gas film _ psia 
P Dimensionless gas film pressure = = 
a 
P. Dimensionless gas film pressure at €= 0, 
see equation (4) 
Pi First order term in pertubation solution, 
see equation (3) 
oa Average value of Po see equation (31) 
P. Ambient pressure psia 
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0 
rosy 


Gas film pressure at [= 0 


Supply pressure 


Slope of Po’, see equation (9) 


Constant defined by equation (35) 
Bearing radius 


Gas constant 


Absolute temperature 
P 
Pressure ratio = = 
a 
Total bearing load 


Circumferential and axial coordinates 
for gas film 


Variable in 2nd part only, defined by 
equation (35a) 


Orifice coefficient 


Defined by equation (17), lst part, and by 
equation (43) - (44), 2nd part 


Bearing eccentricity ratio 


Dimensionless circumferential and axial 
coordinates for gas film 


Constant defined by equation (34) 
Constants defined hy equation (22) - (23) 


Constants defined hy equation (57) - (58) 


2 


Compressibility number = Stas 


() 
Feeding parameter = oo y Bz. 


Gas viscosity 
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peia 


psia 


lbs: sec 


SAAR te a 


a a 


& -<o un 


Bearing length-to-diameter ratio = _ 


Attitude angle 
Angular speed of journal 
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rad/ sec 


APPENDIX 


The Combined Hydrostatic Hydrodyname Gas Journal Bearing 
An Approximate 1st Order Pertubation Solution 

The journal bearing is pressure-fed from orifice restricted feeding holes 
around the circumference in the centerplane of the bearing. For a sufficiently 
large number of feeding holes these may be approximated by a line source 
such that the gas feeding becomes a boundary condition to the differential 
equation describing the flow in the bearing (Ref. 1). The differential equation 
is Reynolds equation which for a perfect gas under isothermal conditions may 
be written: 


& [pH] +8 [pi $8] = uo 


To make dimensionless set: 


ort oe pf h=be I+¢-Cose 
© BlpP $8] +hlpw ae) = ASP 


mA eB 


A let order pertubation solution with respect to the eccentricity ratic € around 
€ 2 0 shall be tried. Therefore set: 


(3) P= R+eP te*R+--- 
and 
P*: B+ 2e RP + --- 
Due to symmetry P. is a function of U only, i.e., e 20 : 


Thus the sero order equation becomes: 
e 


15 


nc ie 


or 


(4) Re I+ p(t) g+ 5 


2 
Since 25 =0 ‘the 1st order equation becomes: 


(5) GH + dH - 4 Bs Ap cing 


(6 H= BP 


The boundary conditions for eq. (5) are: 
a) @*9+27 H(e,t)= H(e+2n,t) 


bp Use H(e,=)=0 


o) CeO mass flow through feeding holes = mass flow through bearing 
For a supply pressure of PP a downstream pressure of P. gas constant R, 
absolute temperature T, orifice radius a, orifice coefficient a and :m orifices 


in series the mass flow M for one feeding hole is approximately: 


M = Tota ig 


For N feeding holes around the circumference the flow per inch is tk 
Equating this flow to the flow through the bearing yields: 


-2 re (3), ~ Ags ae 


Substituting dimensionless parameters gives: 


OH) = AR 


where 


A, a ‘a 


ié 


and 


(8) V= a 


Introducing e4. (3) and (4) yields t and the boundary conditions at [20 


fe AW 


or 


(9) t° Bs [1+ Vis OP I 


and 


(10) Ge). = ¥ cose + #4 Hi 


Approximate solution of differential equation 
To evaluate eq. (5) set 


(11) H =h,(T) sine +h, (0) cose 
which gives: 

(12) 44 “(140 RG *-AB 
where: 

(12a) G= ho) + th,(c) 


Eq. (12) is not readily solved considering that Po is given by eq. (4). Even 


if an exact solution could be obtained it would be very complex and therefore 


be outside the scope of the present analysis which is approximate anyway. 


In the first part 4 


Instead the analysis will be performed in two parts. 


is set equal to a constant. In the second part this assumption will be 
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investigated. 
Rewrite eq. (1): 


$ [8 91 +d(8 81 = Sn + PB] 


which becomes linear in p* when P is a known function. In order to arrive 
at a manageable solution we set P equal to a constant, unspecified for the 
time being. Furthermore, only the lst order pertubation solution is consi- 
dered which is derived by substituting eq. (3): 


dH + 4H -4 qu =-4-P sing 


Proceeding as earlier we get: 


#8 - (sid) G--4e 


with the particular solution: 


(13) G et aan (\-i#) 


The solution to the homogenous equation is: 


(14) G= gull) +i gals) 


where: 
(15) g,(T)=  (arcoshal +csinhal) cospl —(dcoshal +b sinhal) sinft 
(16) gall) = ( beoshal +d sinal) cos BE + (ccoshal + 4 sinhal) sinbf 


a7 gts helt ae =i 


(see Ref. 1) 


i 


Thus the complete solution is: 
(18) Hep a (sine - # cose) +9,(B)sine +9,(6) cose 


Boundary Condition at [=e&: He 
From the boundary condition we get: 


ail8)= ate 


9, (6)= 1+( 
which gives: 
C2, -¥,a~%b 
(19) 
d= d, +ma-ax,b 
where 


- 2h F coshat-sinBt —sinba conte 
2) Com To ( Cosh 2ak~cos2 AF 


: 24 DF soba espe + cosh: sings 
(21) do 1+( Cosh 2u§ ~cos2A§ 


2 —Sinh2at 
(22) Cosheu§—cosepg 
» —sin2Ag__ 
ae % Cosh 2uk = cos 2pg 


Boundary Condition at [20 , eq. (10 
Introducing eq. (18) into eq. (10) and setting [=O gives: 


2 2 
(24) a= +d%+ TA tAG-Bdal +1 Bx- *+Peatdds— 
[4 +ax-+ fa) +[px,-ax,) 
is +Botad.— #)—[pae-aa ]-a+0c 8d] 
(25) ba 
[44 tan,+ Bul +Tpx.-ax,) 
where 


(26) = rae (ER (ip) + gl 


Together with eq. (19) this defines the arbitrary constants in eq. (15) and 
(16), thus completing the solution of Reynolds equation. 


Force Calculation 
The radial and tangential force components are: 


F=-2eR ut Peose do df = ~2R Re " F cose de df 


° 
F.= 2nat{( Psine do dt eRe (t , sing de dl 


where H is given by eq. (18) and Po by eq. (4). Then: 


a1) WERE” Sh tee “l - a( nied 
2a) FELT ap sf Ac 
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The integrals are evaluated by numerical integration. 
The attitude angle i] is given by: 


(29) Q= Arctan (8) 


The total bearing load is given by: 


(30) w= VEtR 


Curves of the attitude angle ) and the dimensionless bearing load mite 
are given in Fig. 1-24forP =A (i.e. P= P,) and for P = Rw. 
where P is the average value of oe and consequently of P. The average 


value of Po is: 


(31) Ray = 


No Feeding, A, =0 


In this case we get: 


~ MiCa-Mede , A . Asinhak-sinB§ —coshat:cospt 
Oe eR THE 


sinhag + cos’ BE 


= MCat Hide. A Accothat cosQt + sinhak- sing 
b +R +2 sinhat + cos'BE 


and therefore c =d=0. Comparing these results to Ref. 1 it is seen that 
a*Tee B and b* A . Thue the preseni analysis reduces to 
Ref. 1 (the lst order pertubation solution) for the case of <A,=0 


No Rotation, A=0 


In this case az2c=z0 and 
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b= — soaks Be stant 
d= — Sahi ,. ca tty 
sinhE Coshé ry sinh& 


Therefore the tangential force F 20, i.e., Q-0 , and the bearing load 


becomes equal to Fe We get: 


= 34¢ sinh(E~ 
aocke ETcoang + Gh sinhé] (es a 


or 


(32) woe * Hrcaa sama (eC vtpe) ipl eh Lporgre) 96] 


where 
Noe * ty 
le) = (eat Gle)=erdle)=  fe*at 
C) ) 
This may be compared to Ref. 3 by means of the following equations: 
Ee 
=f Silt Dae 

U VEO t 


Then eq. (32) may be rewritten: 


sy 3 { <*Cuntse)—wiull 7 Lelia) ecull_} 


TOLRE | EcoshE + ¢ E sinhé 


which is the same as eq. (19A), page 259, of Ref. 3. 


pS 
A Check on the Assumption P = constant 
Since Rr Vit ‘(E-T) is reasonably linear we shall substitute linear 
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A 
expressions for Fe and AP, into eq. (5) and thereby obtain an indication 


of the validity of the previous assumption P = constant. 


Homogenous Equation 
In eq. (12) we introduce the approximation 


(33) ¢ = Al wae +e(I- pene )o] = Aut } 


where 


(34) x= ina 
(35) q° re 


Furthermore set: 


A+0 
(35a) zel+it = | +¢( A+ $) Ayt0 


Then the homogenous part of eq. (12) becomes: 
z 

(36) S +¢2G =0 

with the solution: 


G=Vz [A I,($q2*) +B-1.4(392*)] 


Since Z is a complex variable this expression is rather cumbersome to 


evaluate. Instead the following approximation shall be used: 
(37) G= ZL A-cos(#¢2%) +B-sin($g2*) | 


To evaluate this approximation insert eq. (37) into eq. (36) to get: 


(38) (Uz*-gz)G + gzGu0 
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Comparing the real and the imaginary components: 
I-(Awt 
(39) real part: q= q[I- z ULl+(das ] 


(40) imaginary part: (ax+d)= (Ax+4)[1+ + Serre 


Thus for eq. (37) to be a good approximation the remainders inside the above 
brackets must be small compared to unity. The maximum values of the 


remainders are: 


2 
(41) Max. remainder for real part: ~ Foie rey 


(42) Max. remainder for imaginary part: 77 e ae] 


These equations have been plotted in Fig. 13-20 for a few representative 
values of A, iV and A and thus indicates the regions where the results 
of the analysis is valid. 
Returning to eq. (37) we set: 
$qz" satis 
and 
zs yti b 


From eq. (35a) we get: 


we i [che 
et pre Ag [ (3(K) - + Vie ci? | 
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oe ie) 


(45) x o L(4F) 
7 iC ar 
(46) é= ati (Ay 
where 
(47) + tAy+ i 


Introducing the above equations into eq. (37), setting A = atib and B = ctid 
and using eq. (12a) yields: 


(48) g(t) = al. y Cosa CoshB + d sin sinhB] — b [ d cosa coshp - y Sina sinhp] 
+cly sina-cothB -dcose sinhB] —d[d'sina coshs +y cost sinh] 


(49) 4a (0) = AL d cosa cochB —y sina sinhp] +bLy cosa coshB +d sina sinhB] 


+c[dsina coshB+ycos« sinhAl+d Cy sina cosh ~dcosa siohp] 


Complete Equation 


A particular solution to the complete inhomogenous eq. (12) can be found 
either by using the substitution given by eq. (35a) to obtain a series solution 
or by proceeding with the homogenous solution by means of variation of 
parameters. Neither method yields very practical results. Instead we 
return to eq. (5) and make the following approximations: 
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oe eaten aenee Ba 


i _ £ 
Ro Weak +pli-vie gels Ty-t 
AP= A, -sinh(E-0)+2 sinht 
. 2 SinhE 


Note that for 3 <1 these two approximations for E are nearly identical. 
Then a particular solution of eq. (5) is: 


A 
(50) He Ga Reagan Ae (cheat sinht)| 


Combining eq. (11), eq. (48), eq. (49), and eq. (50) gived the complete 
solution to eq. (5). 


Boundary Condition at [=€ : H=(0 
The boundary condition becomes: 


(51) g(f)= aie rt = aly.com cosh +4, sing, sinhp, ] -bId, cosa, coshf.-y,sing sinh] 
+cLy, sina, Cashfor-d cosa, Sinhp] dd, sine, ceshy+y casa, SinhA] 


(52) ge(€)= = = ALScoset,cosh(,~y, Sina, Sinhp,] + bLy,cosd,cosh/>rtdi sina, sinhf] 
+14, sina ceshB,+y, Cosa, sinh®,]+dLy, sina, cosh®,-d)cose, sinhfy] 
where 4,, P,, y and 6, are derived from eq. (43)-(46) by setting 


tzk (i.e., r* 1). 
Solving eq. (51) and (52) with respect to c and d yields: 


(53) C=Co-A,A=)dQb 


(54) d=d,+),a-),b 
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where: 


_ (55) Co™ (+5; 4.(E)) sing, cosh. + (y,-9,(6)-dyq (8) st 


(y3 +d?)( coshep, —cos2u,) 


2 (€) -d,-a,l€)) sing, coshfy —(x,:q,(8)-+d,90(8) in 


(56) d= 
_ (y?+0;') (cosh2p, -cos2a,) 


= Sin2a, 
(57) A= CosheB,-cosda, 


) d2= Cosh2f,- cos2e, 


Boundary Condition at [20 , see eq. 10) 
Introduce the following constant: 


(59) do = (a), 
(60) BL * (Bless 
61) Je = eee 
(62) &* (dh. 
2) ae (fhe $ fe Keele 


ad 


al=(H6)..0 $e elect) 


(64) 
cs) gh (A =a PLAC pe) He Vile 1/;, (ree 


27 


eos tM aie f 


(66) 


(67) 


(68) 


(69) 


(70) 


(71) 


(72) 


(73) 


(74) 


(75) 


(76) 


(77) 


(78) 


(79) 


&= (8) 24% [s(t psn) Viet 400 ce) /a (tee 


VE Yo Cod. Cosh. +d, sing, sinh 

g= d, cose, CoshAy —y, sind, sinhpy 

3 ™ y, Sina, Coshfo —d Cos, SinhBe 

T= A Sind, Coshp, t yo Cosa, SinhB 

V'® Yo Cosa, CoshBe +dy SindsSinh Be dod +A,T 
g'* d, com, CoshBey, sina, Sinh By — aT — fed 
3'=¥ Sind, Coshe-do cosd, Sinh Py +d,V—fs g 
T'= 5, sing, Cosh ety, CoSd, SinhBs tdeg + ALY 
s= axs 

Y= v-As-d,T 

M*3+)6-),T 

Wev'-do'-a7 


m'= g'+ho' -), t' 
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(80) Q* SBM Ag)+A4] +(6'- Hal -(t'- fra, 
(31) M25 TA [io Ag hs +Aq cosh] +(t’- ferjc,+(6'- Kola, 


Then we get from the boundary condition eq. (10): 


(82) as 


2 


J-o[ —) 


eames (Ay-y7+[B 9-4] 


Force Calculation 


The radial and tangential force components are: 


aw 
(84) FF = -28Re| | # cose do dl 


ry 
(85) F = zane | # sino de dl 


where: 
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(a6) H ={ Gatttag (8-Aq) coset) ~arcoshl) + 9\(0) } sino 
+{ staat (T-Ag)(sinh(-0)+ esinht) + Ga(d)} cose 
(87) Be , I+ 92(§-0) 


and g,(5) and 926) are given by eq. (48) and (49). Then: 


5 
a ; (T-AgXsinh(§-() +a sinh) +44 ( 
(88) For KE" +\; +4E-T) dt 


wn fics $ an a ie aab dt 


Vis PET) 
These integrals are evaluated numerically. The total force is: 
(90) W=VR+R 


and the attitude angle is determined by: 


(91) Q= Arctan (ft) 
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